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We present a new approach to quintessential inflation, in which both dark energy and inflation 
are explained by the evolution of a single scalar field. We start from a simple scalar potential with 
both oscillatory and exponential behaviorWe employ the conventional reheating mechanism of new 
inflation, in which the scalar decays to light fermions with a decay width that is proportional to the 
scalar mass. Because our scalar mass is proportional to the Hubble rate, this gives adequate reheating 
at early times while shutting off at late times to preserve quintessence and satisfy nucleosynthesis 
constraints. 

We discuss a simple model which solves the horizon, flatness, and "why now" problems. Without 
any additional tuning of parameters, this model satisfies all constraints from CMB, large scale 
structure, and nucleosynthesis. The predictions for the inflationary spectral indices are ns = nr = f . 
In this model we are currently beginning the third cosmic epoch of accelerated expansion. 



PACS numbers: 98.80.Cq, 95.35.+d, 98.70.Vc 

If we assume the correctness of the standard Fried- 
mann equation evolution, the existence of dark energy 
engenders two profound dilemmas. The first is the cos- 
mological constant problem, and the second is the "why 
now?" problem. Quintessence models attempt to address 
the second problem by introducing a very weakly coupled 
scalar field whose potential and/or kinetic function have 
special properties. One of the most successful approaches 
to quintessence 0,0] is to combine tracking with an os- 
cillating behavior in the quintessence potential. In such 
models the equation of state parameter w(z) has a peri- 
odic component, leading to occasional periods of acceler- 
ated expansion during epochs where w(z) ~ — 1. 

It is natural in this context to ask whether the 
quintessence scalar could replace the inflaton. The idea 
of quintessential inflation has been examined by a num- 
ber of authors |2|-[lJ|- The straightforward approach is 
to cobble together a scalar potential which has both a 
flat, large vev portion (for inflation) and a flat, small 
vev portion (for quintessence). These features are con- 
nected by a steep step which corresponds to a period of 
cosmic kination. As discussed in |3, LXQl , such models suf- 
fer from generic problems. First, they require significant 
ad hoc tuning to simultaneously produce the features of 
inflation and quintessence. Second, they require a "ster- 
ile" inflaton, in order to avoid the decay of the putative 
quintessence scalar at the end of inflation. This in turn 
requires new mechanisms for reheating, such as gravita- 
tional de Sitter phase particle production, leading to dif- 
ficulties in satisfying the constraints of CMB anisotropics 
and of primordial nucleosynthesis. 

Our approach to quintessential inflation is to take ad- 
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FIG. I: The upper (lower) panel shows the scalar potential 
V(6) (V(<j>)) for 6=1, plotted on a logarithmic scale. 



vantage of the tracking and oscillatory potential features 
that work so well in addressing the "why now" problem 
of quintessence alone. We will describe a model with 
a simple scalar lagrangian with exponential and oscilla- 
tory features. The model uses the conventional reheat- 
ing mechanism of new inflation , in which the scalar 
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decays to light fermions with a decay width that is pro- 
portional to the scalar mass. We show that the scalar 
mass is proportional to the Hubble rate. As a result, 
the model has adequate reheating at early times while 
naturally shutting off at later times. 

We present a successful model with only three ad- 
justable parameters. One parameter controls the period 
between inflationary epochs, a second controls the over- 
all decay width, and the third parametrizes our ignorance 
about the relative fraction of matter versus radiation pro- 
duced by reheating. These three parameters are adjusted 
to produce sufficient inflation along with the correct frac- 
tions Q r /Cl\, £l m /Q\ of radiation, matter, and dark en- 
ergy, as measured today. Having thus fixed the model we 
find that we automatically satisfy all constraints of pri- 
mordial nucleosynthesis, CMB, and large-scale structure. 



AN OSCILLATORY POTENTIAL 

We start with a simple model with a single real scalar 
quintessence field 9. The action is 



d X\J—g 



f{9)g^d^9d u 9 - V(6) 



(1) 



where the kinetic function f(9) and potential V(9) are 
given by: 



m = -^sin^; 
V(6) = pocos 2 0exp 



- (26> - sin 26>) 
b 



(2) 
(3) 



where Mp is the Planck mass: 1.22 xlO 19 GeV; po is 
the dark energy density observed today: ~ (10 -4 eV) 4 ; b 
is a dimensionless parameter which controls the periodic 
behavior. A canonical kinetic term can be restored via a 
field redefinition 9{x) — > 4>(x), where 



4>{x) — 0o cos# 



(4) 



with 0o = y/3M£/nb 2 . 

In the approximation where we ignore the energy den- 
sity of radiation and matter, and where the only friction 
is from the metric expansion, the evolution of the model 
can be solved analytically. Energy conservation requires: 



-3H(l + w)p e 



(5) 



where H is the Hubble rate, w is the equation of state 
parameter, and pg is the dark energy density. The solu- 
tion to this equation as a function of the scale factor a(t) 
is: 



pe(a) = p exp 



„ f a da , ... 
-3 / — (l + w(a)) 
Ji a 



(6) 



We also know that 
V{6) 

Making the Ansatz 



1 



(1 - w)pe 



w(a) = —cos29(a) 



(7) 



(8) 



one immediately gets a solution to the (flat) Friedmann 
equation combined with the relations (|6I7|I : 



Ha) = lna ; 
)(a) = —cos [b In a] 



(9) 
(10) 



The expectation value of the quintessence field 9 
evolves logarithmically with scale factor from a positive 
initial value to zero today. The potential V(9) has the 
"Slinky" form 1 shown in Figure^ Accelerated expansion 
corresponds to epochs (such as today) where 9 is evolv- 
ing through one of the flat "steps" of the potential. From 
(|2*|) we see that the kinetic function is simultaneously sup- 
pressed in this epochs, slowing the roll of the scalar field 
evolution. Note also that our potential V(9) — ► as 
9 — > — oo, which corresponds to t — > oo; this is as desired 
for a quintessence model. 

We can understand the same behavior by looking at 
the "canonical" scalar <fi in (0J. The potential V(<p) re- 
sembles a series of descending ramps, with the value of 4> 
bounded between — (f>o and <f>Q. Shifted around a generic 
vev ((f)), the scalar has an effective mass-squared given 
by 



H 2 



4 



cos 3 9 



sine 



- 2sin 2 26> 



-sin 2 26>, (11) 



where H is the Hubble rate. Instead of a suppression of 
the kinetic function when 9 is a multiple of tt, the scalar 
(j) has an effective mass that diverges when ip = ±0o, i.e., 
at each junction of the descending ramps. 

The equation of state parameter w(a) has the same 
periodic form assumed in two recent phenomenological 
analyses The analysis in 01 showed that, for 

6=1, w(a) as given in l|10|l is consistent with all current 
data from observations of the CMB, Type IA supernovae, 
and large scale structure. It follows a fortiori that our 
model with any choice of b less than one also agrees with 
this data. 

To complete the model, we will assume that the 
quintessence field cf> has a weak perturbative coupling to 
light fermions. This is the standard reheating mechanism 
of new inflation [ll|. To avoid the strong constraints on 



where po is the dark energy density at a = 1 (today). 



1 A Slinky™ is a large spring with a very weak spring constant. A 
Slinky will "walk" down a staircase, producing in time-lapse a 
pattern resembling the top panel of Figure HI 
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long-range forces mediated by quintessence scalars [Tij . 
it is simplest to imagine that our scalar only has a direct 
coupling to a sterile neutrino. This is sufficient to hide 
the quintessence force from Standard Model nonsinglet 
particles , while still allowing the generation of a ra- 
diative thermal bath of Standard Model particles from 
quintessence decay. 

With this assumption for reheating the evolution equa- 
tions for quintessence, radiation, and matter become: 

pe = -3H(l + w)p e - k m^(l + w)p e ; 

p r = -4ff(l + w)p r + (l-/ m )fc m^(l + w)p g ;(12) 

Prn = -3iJ(l + Vj)p m + fmkom^l + w)p 6 \ 

where fco and f m are small dimensionless constants. From 
we see that, as long as 9 is not near a multiple of 7r/2, 
it is a reasonable approximation to make the replacement 



k m 4 



kH 



(13) 



where k is another small dimensionless parameter. This 
replacement decouples the 9 evolution equation from the 
Fricdmann equation, giving an immediate analytic solu- 
tion: 



pe(a) = p exp 



k) (29 -sin 26») 



(14) 



where 0(a) and w(a) are still given by lj§)l- illO|) . 

We have used this approximation in the solutions 
quoted below. Clearly this approximation overestimates 
the reheating effect for epochs where 9 is close to an odd 
multiple of 7r/2, however this is a small effect since these 
are the epochs of maximum radiation or matter domina- 
tion. The approximation also breaks down in the epochs 
of maximum inflation, i.e. when sin9 — > 0. However 
when sin 8^0 we also expect strong-coupling physics 
to enter in a full theory, providing a cutoff for the naive 
vanishing of the kinetic function of 9. Since we cannot 
compute this effect without a full Planck-scale theory, we 
may as well stay with the approximation (|13fl . 



RESULTS 

Figure and Table 1 show the results obtained from 
our model with b = 1/7, k = 0.06, and f m = IO -11 . 
Shown are the relative energy density fractions in dark 
energy, radiation, and matter, as a function of log a. We 
have chosen to integrate the evolution equations starting 
from a — 10~ 42 , which in our model corresponds to a 
temperature of slightly less than 10 16 GeV, and an initial 
comoving Hubble radius of about 100 Planck lengths. For 
simplicity we have also chosen the value of w(a) now to 
be exactly —1. Neither of these choices corresponds to a 
necessary tuning. 

Our three adjustable parameters have been chosen 
such that the values of Q r /£l\, f2 m /f2A come out to their 
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FIG. 2: History of a slinky universe. Shown are the relative 
energy density fractions in radiation (green), matter (blue), 
and dark energy (red), as a function of the logarithm of the 
scale factor a(t). 
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TABLE I: The relative density fractions of dark energy, ra- 
diation, and matter, as a function of the scale factor. The 
future projections are only valid if the quintessence scalar is 
still decaying. 



measured values at a = 1, and such that we have suffi- 
cient inflation. The latter is checked by computing the 
ratio of the fully inflated size of the intitial comoving 
Hubble radius to the current comoving Hubble radius. 
This ratio is about 3 in our model, indicating that the 
total amount of inflation is indeed enough to solve the 
horizon problem. The flatness problem is solved because 
the total Q r + Q m + J7a = 1 within errors. 

From Figure (J2J we see that we are currently beginning 
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the third epoch of accelerated expansion. The first epoch 
of inflation accumulated about 18 e-foldings. Quantum 
fluctuations during this epoch produced the spatial in- 
homogeneities responsible for large scale structure and 
CMB anisotropics observed today. Constraints on the 
physics responsible for the primodial power spectrum of 
these density fluctuations can be set with WMAP and 
2dF data, under the assumption that the Hubble rate is 
dominated by the contribution from during the ob- 
servable part of inflation 0] . The scalar and tensor per- 
turbations with comoving wavenumber k are solutions of 
the equation 



small current mass of ~ 10~ 36 eV. Even with such an 
effect, our model satisfies the constraints of successful 
nucleosynthesis [2cj . 

Figure J2J) also shows one (very brief) prior epoch of 
matter domination. This only occurs if there is some 
suitably heavy and long-lived matter around to go out of 
thermal equilibrium when a ~ 10 -22 , e.g. a superheavy 
neutrino. 



REMARKS 



< + 



k z )u k = , 

z 



(15) 



where prime denotes derivative with respect to conformal 
time t, and z — acj>/H for scalar perturbations while 
z = a for the tensor case. The power spectra have the 
form 



Vn(k) 



where e 



H 2 



H 2 f k . 



fk 



k 3 uj 
a 2 H 2 



(16) 



-H/H 



For standard inflation models 



(l?|. the solutions of ifTKjl have the property that 
a 2 H 2 /k 3 for small wavenumbers k/aH — > 0, and thus 
fk — ► constant. This gives the standard formulae for the 
spectral indices ns — 1 = — 4e — 26 and rpp — 1 = — 2e, 
where 5 is defined by 5 + e = e'/ (2aHe). 

In our model (|15|l can be solved in the limit k /aH — > 
by changing variables from r to 9: 



d2u k _ 2(1 _ o ej n 2 m duk 

d02 b {i osm o) de 
c(9) = fcot6» + 1 - 

m = 



c{6)u k = 



6 (cos 26- 1) 



8 , 12 



§sin26> 



sin 



(17) 
(18) 
(19) 



where the first/second expression for c(9) is for the 
scalar/tensor case. The solutions which match to the 
correct large k behavior are: 



N a 

—=— sm( 

Vkk 



N a 



(20) 



for the scalar and tensor cases respectively, where N is a 
constant. 

We can plug these solutions into the expressions i|16|) 
for the power spectra, using also that, for our model, 
e = 3 sin 2 6*. Remarkably, one finds that ns — ut = 1. 

To compare with contraints from primordial nucleosyn- 
thesis, we note that the second epoch of dark energy dom- 
ination ended well before a ~ 10~ 10 , the time at which 
nucleosynthesis occured. Indeed dark energy reheating 
effects are completely negligible from a ~ 10~ 10 until 
today. A possible exception is a small reheating effect 
occuring at very late times in our current epoch of ac- 
celerated expansion; this only occurs if the quintessence 
scalar is still managing to decay despite its incredibly 



It is fair to say that, compared to the simple model pre- 
sented above, the standard new inflation scenario looks 
rather extreme. In the evolution history portrayed in Fig- 
ure [2 the interplay between inflationary expansion and 
reheating is much milder. In fact, apart from a few very 
brief periods, reheating effects in our model do not ac- 
tually increase the temperature of the thermal radiation 
bath. Instead, the temperature is almost always decreas- 
ing, but it decreases more slowly than in the standard 
ACDM evolution. 

Our particular realization of Slinky inflation must be 
regarded as a toy model (pun intended), since <pQ > Mp 
and we have not invoked any consistent Planck scale 
framework such as string theory. However the model it- 
self is surprisingly simple, and the physical picture which 
emerges from it has some compelling features. These are 
worthy of further investigation. Also, since the target 
space parametrized by our scalar is S 1 , it would be in- 
teresting to extend this scenario to a class of nonlinear 
sigma models with other compact target spaces. 
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